Escape configuration lattice near the nematic-isotropic transition: Tilt analogue of 

blue phases 
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We predict the possible existence of a new phase of hquid crystals near the nematic-isotropic {NI) 
transition. This phase is an achiral, tilt-analogue of the blue phase and is composed of a lattice of 
double-tilt, escape-configuration cylinders. We discuss the structure and the stability of this phase 
and provide an estimate of the lattice parameter. 
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Liquid crystals are soft materials; this makes the exis- 
tence of phases such as the blue phases 0, 01 , the twist 
rain boundary phases 0, 0| and the smectic blue phases 
'1] possible. In this paper we explore the possible 
existence of a new phase of liquid crystals intervening 
the nematic {N) and the isotropic (/) phases. The pro- 
posed phase is composed of escape-configuration cylin- 
ders Slid can be thought of as an achiral, tilt ana- 
logue of blue phases (see FigQ. Blue phases are com- 
posed of double-twist cylinders and intervene between 
the cholesteric and the isotropic phases. Interestingly, 
the stability of double-twist cylinders forming blue phases 
is not solely due to the chirality of molecules. Double- 
twist cylinders are favored over the simple cholesteric 
with twist along a single axis if the coefficient of the 
saddle-splay deformation term in the elastic free energy 
is positive [llQ. The saddle-splay deformation (see (0)) 
corresponds to a surface-like term in the free energy (it is 
a total divergence and integrates to the bounding surface 
by Gauss's theorem). Therefore blue phases necessarily 
require the presence of interfaces for their stability. The 
gaps between the double twist cylinders in the cubic blue 
phases are believed to be filled with higher symmetry, 
lower order isotropic material^]. 

For a two-dimensional nematic confined within a cir- 
cle, and subject to the boundary condition that the ne- 
matic director is radial at the circumference, it is a topo- 
logical imperative that the director field inside the cir- 
cle have point singularities (disclinations) with a total 
charge of +1. However, a similar topological require- 
ment does not apply to a three-dimensional nematic con- 
fined within a cylindrical capillary with radial bound- 
ary conditions. The -1-1 line singularity can be removed 
by the escape of the nematic director into the third di- 
mension The singular disclination line is stable if 
ln(i?/a) < f (ifa/isri)^/^^ where R is the radius of the 
capillary, and a is a microscopic cut-off (see ^ below). 
Clearly, escape-configurations are favored over configu- 
rations with topological singularities for reasonable cap- 
illary radii unless the splay elastic constant is unusually 
small. The director field n in the escape-configuration 
cylinders (FigEJ is tilted with respect to the axis of the 



FIG. 1: Schematic of the antiferroelectric square lattice of 
double-tilt cylinders. As in the case of the cubic blue phases, 
the gap between the cylinders is filled with the higher symme- 
try, lower order isotropic phase. The director field maintains 
continuity across the planes joining the axes of the neighbor- 
ing cylinders. See Fig.Q for a two- dimensional representa- 
tion of a single double-tilt cylinder. 



cylinder, and has splay as well as bend deformation. This 
tilt of the director can be resolved in two mutually or- 
thogonal directions in a plane perpendicular to the axis of 
the cylinder; they are double-tilt cylinders. These cylin- 
ders clearly have a non-zero saddle-splay deformation: 
the two principal curvatures of the nematic director field 
at any point in the cylinder have opposite signs. 

Can a phase composed of the double-tilt cylinders ex- 
ist near the nematic-isotropic (NI) transition for reason- 
able values of parameters? In what follows we answer this 
question in the affirmative, and propose a possible candi- 
date for such a phase. In the spirit of the "low chirality" 
analysis of blue phases P, ? we use the continuum de- 
scription of nematics and ignore the effects of thermal 
fluctuations. As is the case with the low chirality theory, 
the theory presented in this letter is not a theory of phase 
transitions based upon an order parameter description. 

Our results are summarized in Figs.l^ 21and|SJ|. The 
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FIG. 2: Schematic of the escape configuration. The solid 
vertical lines represent the nematic-isotropic interface. The 
angle 9 preferred by the weak anchoring energy differs from 
the angle (j) which minimizes the total energy. The nematic 
director is tangential to the curved lines. 



lattice depicted schematically in Fig.lQ]) is stabilized by 
the competition between (i) the saddle-splay deforma- 
tion ||2J); and (m) the bulk nematic deformation the 
nematic-isotropic interface energy |(2Jl, and the weak- 
anchoring energy at the interface In the structure 
proposed in this paper, the gaps between the cylinders 
are filled up with the high symmetry, low order isotropic 
material. Notice that the lattice shown in Fig.QJ is an 
antiferroelectric lattice: neighboring cylinders have op- 
posite escape directions to maintain the continuity of the 
director field across the planes joining the cylinder axes. 
Each double-tilt cylinder can be characterized by the di- 
rection of escape (up or down), and is a polar object 
even if the nematogenic molecules are not polar. Fig.Q 
shows that the proposed phase has a small region of sta- 
bility near the NI transition temperature in the space 
of the dimensionless parameters ^24 — K24,/ K, the rel- 
ative strength of the coefficient of the saddle-splay term 
to the typical coefficient of the bulk Frank free energy, 
and A: = 7/(7, the relative strength of the coefficient of 
the nematic-isotropic interfacial tension to the coefficient 
of anchoring energy at the interface. In the case of cubic 
blue phases, the pitch of the cholesteric provides a natu- 
ral length scale, and the nematic director twists through 
7r/4 from the center of the double- twist cylinder to its 
boundary. This ensures that the director is continuous 
across neighboring double-twist cylinders. In the present 
problem, there is no "natural" length scale at the lowest 
order of the continuum description that we employ. We 
therefore calculate the radius of double-tilt cylinders by 
solving the free boundary condition that arises naturally 
from the variational problem of minimizing the total free 
energy (see Fig.lO shows the variation of the radius 
of the cylinder (which is half the lattice parameter) with 
/i24 for various values of k. 

In order to derive these results we begin by listing the 
contributions to the total free energy. The bulk Frank 



free energy density is 

/F = ^(V.n)^ 



-(ri-(V xn)f 



K3, 



^(nx(Vxn))2, (1) 

where Ki,K2j and K3 are respectively the splay, twist, 
and bend elastic constants. In what follows, we use the 
"one-constant" approximation Ki = K2 = K3 = K, 
which is indeed a reasonable approximation near the 
nematic-isotropic transition temperature j^. 

In addition to the bulk contribution, we have the 
surface-like saddle-splay elastic free energy density 

/ss = i^24 V • [n( V • n) - (n • V)n] . (2) 

We note that the sign of K24 is not dictated by stability 
conditions. 

The free energy associated with the nematic-isotropic 
interface is composed of two parts, the usual surface ten- 
sion part 
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(3) 



and the Rapini-Papoular [loj "weak anchoring energy 



(4) 



where 6 is the preferred angle of the nematic director at 
the interface (see Figl^J, and the integrals are over the 
interface. Experimentall y, v alues of 9 as large as about 
45° have been reported The spontaneous nonzero 

tilt of the director at the interface has its origin in 
the order electric polarization |l3| arising from the gra- 
dient of the nematic order parameter across the nematic- 
isotropic interface. Since the present model is not based 
upon an order parameter description, it is sufficient for 
our purpose to recognize that 9 is nonzero and use @ 
for the energy cost of deviations from the angle 9 at the 
interface. 

We now discuss the energetics of a single escape- 
configuration cylinder at the NI transition (at the transi- 
tion point, there is no free energy difference between the 
nematic in the escape cylinders and the isotropic material 
filling in the gaps between them). The total free energy 
for our model is 

Ft — Fp + Fss + FsT + Frp, (5) 

where Fp = J fpd^x, and Fss = / hs<Px. 

We parameterize the nematic director as rir = 
sin ri(r), = Q,nz — cosr2(r) in cylindrical polar co- 
ordinates, where ^{r) is the angle made by the director 
with the axis of the cylinder (see FigEl). The free energy 
per unit length of the cylinder of radius R is 
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FIG. 3: Variation of the free energy per unit length of a 
single double-tilt cylinder as a function of 4> for 9 = 7r/4, 
k — 7/cr = 1; /i24 — —1.3 (dashed line), /i24 ~ —1.373 (dash- 
dotted line) and /i24 = —1.5 (dotted line). For /i24 — —1.5, 
/t has a minimum. 



where K = + and the length along the 

cylinder axis is unity. The Euler-Lagrange equation reads 



IdVL 1 

— ; 7T smil COS 52 = 0. 



(7) 



dr'^ ' r dr 

The variational problem of minimizing the total free en- 
ergy leads to the free boundary condition 

dVL^ k , o- , M 

K — \r=B. H sin i7 COS iZL^fl H — sin2(S2 — 9)\r=R ~ 0, 

dr r 2 

(8) 

where R is the radius of the escape cylinder. This con- 
dition which ensures torque balance at the boundary of 
the escape cylinders picks out the radius of the escape 
cylinders and hence the lattice parameter. 

We recognize that the angle (f) made by the nematic 
director at the NI interface is not necessarily equal to 
the angle 9 preferred by the anchoring energy. Therefore, 
we first solve Eq.© subject to the condition that the 
nematic director makes an angle (j) at the NI interface, 
and then proceed to minimize the total free energy with 
respect to (j) (see Fig. 

The solution to the Euler-Lagrange equation is 0| 
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The free boundary condition relates the radius R 
of the cylinder to the angle Q. — (f) aX the boundary of the 
cylinder, so that 
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(6* - 0)(2sin(/)-f ^sin20), (10) 



where A = K/a is a length, fi — 2/i24 + 1, with /i24 — 
K24/K. Substituting the expression ffT^ for the radius 
in /t yields the energy per unit length (measured in units 
of K) of a single double-tilt cylinder of radius R: 



/t = 2^ 



. 2 ^ 
sm — 
2 



fj, sm 



(11) 



FIG. 4: "Phase diagram" in the k — fi24- plane for 9 = tt/A. 
Escape cylinders with a finite radius are stable in the region 
labelled I, and unstable in regions labelled II. The range of 
stability is larger for smaller values of k, i.e. for stronger 
anchoring and weaker nematic-isotropic surface tension. The 
continuum description holds in the hashed region (see the 
discussion following Eg. lllU 'l. 
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FIG. 5: Radius of the cylinder R measured in units of A24 = 
K2i/(T (which is about 0.3 ^m for a typical nematic) as a 
function of /^24 for 9 = n/A; k = 0.5 (circles) and k = 1 
(triangles). 



where A; = 7/0-, and we have defined f{9, </>) — esc 2(0 — 
(j)) (sin 4> + 11 sin 2(f)) . We numerically minimize the free 
energy (|ll|l with respect to cj) to find the dependence of 
R on /i24 (Fig.®), and of on (Fig.®). 

We now estimate the approximate value of the radius 
at which double-tilt cylinders are energetically favored 
over the uniform, undistorted nematic at the NI tran- 
sition. With fi24 ~ —1.5 and A24 w 3 x 10^''' cm; the 
radius of the double-tilt cylinders is about 0.04 /itm for 
fc w 1, whereas it is about 0.1 /xm for k « 0.5. It is im- 
portant to note that the continuum description that we 
have employed fails at very small radii of the cylinders. 
The NI transition is first order in nature, and the order 
parameter coherence length ^ is of the order of molecular 
dimensions (see 0). The cylinder sizes estimated above 
(for a reasonable range of parameter values) are about 
ten times the order parameter coherence length. From 
Fig.® we expect the continuum description to be valid 
within the shaded region indicated in Fig.®. 
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FIG. 6: The variation of (j) with /i24 for k — 0.5 (circles), 
k = 1.0 (triangles), for 6 — n/4. As expected, (j> is closer to 9 
for smaller values of k and /i24- 



Fig.Q) depicts the range of stability (labelled I) of the 
double-tilt cylinders in the fc — /i24 plane. In the region la- 
belled II on the right- hand side of the spike in Fig.Q the 
gain in energy from the saddle-splay term is not enough 
to overcome the energy costs from the Frank free energy, 
the surface tension and the anchoring energy terms. On 
the left hand side of the stable spike in figure 01 as the 
signed value of K24 decreases, and (p at the boundary ap- 
proaches 9 in order to gain energy from the saddle-splay 
term (see Fig.®). However, the torque balance condi- 
tion ijHll at the boundary forces the relation , thereby 
reducing the radius of the cylinders (see Fig.lISJl). Equa- 
tion pU|) implies that the radius of the cylinders goes to 
zero for 1 + /i cos 0^1. Finally we discuss the possi- 
bility of forming lattices of escape cylinders. For a two 
dimensional lattice, continuity of the director field across 
the planes connecting the axes of neighboring double-tilt 
cylinders enforces the condition that these have oppo- 
site escape directions. Clearly, the triangular lattice is 
frustrated. The proposed square lattice is antiferroelec- 
tric. Having demonstrated the feasibility of forming a 
two-dimensional lattice of double-tilt cylinders for rea- 
sonable parameter values, we note the following: (i) it 
may be possible to construct three dimensional lattices 
out of escape configurations which are not necessarily 
cylindrical in shape. We have explored the possibility 
of forming three dimensional structures of right circular 
double-tilt cylinders. Cubic lattices (BCC and FCC) can 
be formed if the angle (p = tt/A. For all other angles 
we find that matching the director field at the bound- 
aries of neighboring cylinders introduces additional dis- 
tortions in the director field of the double-tilt cylinders. 
Energetics of these three dimensional structures is be- 
yond the scope of this paper, (ii) It is possible to con- 
struct lattices composed of hyperbolic escape configura- 
tions (which correspond to escaped -1 disclinations), and 
{Hi) the double-tilt cylinders can form a molten line liq- 
uid analogous to the molten flux line liquids^] and the 
chiral line liquids 16, of liquid crystals. In conclusion 
we have clearly demonstrated that the structure proposed 



in this paper has a lower free energy than the nematic 
as well as the isotropic phase near the NI transition. 
Whereas we do not claim that this structure is the global 
minimum of the free energy, our demonstration strongly 
suggests the possibility of lattices composed of double- 
twist cylinders. 

Freeze fracture electron microscopy of pure compounds 
exhibiting a broad NI transition may be a possible way 
to search for the proposed phase, since the length scale in- 
volved (0.05/xm to O.I^to) is inaccessible to other exper- 
imental probes. We urge experimentalists to look closely 
for signatures of such structures near the NI transition. 
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